Effective field theories for disordered systems from the logarithmic derivative of the 
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We consider a spinless particle moving in a random potential on a d-dimensional torus. Introduc- 
ing the gradient of the logarithm of the wave-function transforms the time independent Schrodinger 
equation into a stochastic differential equation with the random potential acting as the source. Us- 
ing this as our starting point we write functional integral representations for the disorder averaged 
density of states, the two point correlator of the absolute value of the wave-function as well as the 
conductivity for a d-dimensional system. We use the well studied one dimensional system with 
Gaussian disorder to illustrate that these quantities can be computed reliably in the current for- 
malism by using standard approximation techniques. We also indicate the possibility of including 
magnetic fields. 
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I. INTRODUCTION 

The study of the localization of non-interacting elec- 
trons in disordered media has progressed considerably 
since the initial work of AndersonEl in 1958. Early work 
on one dimensional (ID) disordered systems iiicliided the 
study of the spectral densities- by HalperiroJ'ta as well 
as a calculation by BerezinskiiEl using diagramatic tech- 
niques showing that all states are localized in ID dis- 
ordered systems, although this is generally difficult to 
extend to higher dimensions. Abrahams et al.u intro- 
duced a scaling theory of localization predicting that a 
metal-insulator transition occurs in dimensions greater 
than two, although there seems to be experimental ev- 
idence for a transition in two dimensionso. Making use 
of the replica tricmL the problem was mapped onto a 
non-linear cr-modelBli30, which gave quantitative con- 
firmation of th e spa ring approach. Efetov's supersym- 
metry approachE3E3 introduced a mathematically more 
rigorous alternative to the replica trick, which he used 
to prove, amongst other things, a conjecturCj-of-Gor'kov 
and EliashbergJi-3 that random matrix theorytatij can be 
applied to the energy level statistics of particles in disor- 
dered systems. 

Notwithstanding the considerable amount of work that 
has gone into the investigation of the localization prob- 
lem, there are still many outstanding-problems, for in- 
stance the lack of a order parametenl3 to describe the 
2nd order metal-insulator phase transition. Also, finding 
an analytically tractable description of the localization 
problem, of which there has been little progress, would 
lead to a better understanding of disorderecLbased phe- 
nomena, such as the Quantum Hall Effectta. For this 
reason, any additional aproaches for studying the local- 
ization problem, possibly leading to new insights, are use- 
ful. 

In general, we would like to calculate disordered aver- 
ages of observables that depend on a random potential 
V(x). These disordered averages can be calculated when 
the exact dependence of the observable on the random 
potential is known. However, when this dependence is 



not known, for example, the density of states and corre- 
lators of the wave- function, other methods of averaging 
these observables over the disorder are needed. Usually, 
the disorder averages of advanced or retarded Green's 
functions, G ± (E) = (E - H ± ie)" 1 , are calculated since 
their dependence on V(x) is known. These averages are 
then related to the averages of the observable. Thus, 
one would calculate the average of the advanced Green's 
function and then relate it to the density of states using 

(p(E)) = -l\m^- 1 lmTT(G(E)), (1) 

where the angle brackets denote averaging over the dis- 
order. Both of the main field theoretic techniques for in-, 
vestigating-disordered systems, the supersymmetryPIF 3 ! 
and replicao methods, are based on calculating the aver- 
ages of products of Green's functions using a generating 
function and then extracting information from the result. 

In this paper we would like to propose a complemen- 
tary approach for calculating disorder averages. This ap- 
proach entails a transformation where we change from 
the random potential V(x) to a new set of random vari- 
ables, which can be related to the logarithmic derivative 
of the wave-function and energy of a particle moving in 
the random potential. Using this transformation allows 
us to calculate directly averages of the density of states 
and correlations of the wave-function and its absolute 
value. 

In section |J, we will introduce the formalism, both for 
one-dimensional systems and for higher dimensions, and 
show how disordered averaged observables are calculated 
within this framework. Since the one dimensional sys- 
tem with Gaussian disorder is probably the best studied 
disordered-system, with a variety of well known results 
availableBE3, it is ideal for testing and developing ap- 
proximation techniques within our formalism with the 
ultimate aim of extending these techniques to higher di- 
mensions, and possibly also to the case of a magnetic 
field. Therefore we focus in section III on the one dimen- 



sional Gaussian disordered system in order to illustrate 
how the formalism can be applied, using standard ap- 
proximation techniques, to recover known results for the 
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density of statesa, and to obtain results for the 2-point 
correlator of the absolute walue of the wave functioned 
as well as the conductivity^. In section III we also give 
a realization of the model to show how the parameters 
describing the Gaussian disorder can be related to micro- 
scopic quantities. In the final section, section IV, we nu- 
merically calculate and generate plots of the main results 
obtained in section [II in order to obtain a understanding 
of the results. 



II. FORMALISM 

We consider a particle moving on a d-dimensional 
torus, and in a periodic, random potential caused by 
impurities in a system. We wish to calculate observed 
quantities of this particle, which under the assumption of 
self-averaging implies the averaging of these observables 
over the different realizations of the random potential, 



(6> = z- 1 / [dv]d[v]P[v], 



(2) 



where Z = J[dV]P[V] and P[V] is the probability dis- 
tribution describing the random potential. If we as- 
sume that the impurities are quenched, then the move- 
ment of the particle is described by the time indepen- 
dent Schrodinger equation. We impose periodic bound- 
ary conditions and for the moment assume that time re- 
versal symmetry is not broken, so that the wave function 
can always be chosen real. 

To introduce our formalism, we consider the loga- 
rithmic derivative of the wave functionO instead of the 
wave function, and correspondingly change from the 
Schrodinger equation to the equation of motion of this 
new variable. There are several advantages to this 
transformation, particularly from a functional integra- 
tion point of view. Firstly, in contrast to the Schrodinger 
equation in which the random potential multiplies the 
wave function, the equation of motion governing the new 
variable is a non-linear stochastic differential equation 
in which the random potential simply plays the role of a 
random source. Secondly, the physical irrelevant normal- 
isation of the wave function is eliminated. Lastly, as is 
well knownta, the localization length cannot be extracted 
from the correlations of the wave function as these are 
always short ranged due to the random phase cancella- 
tions. Instead, the correlations of the absolute values of 
the wave function should be computed. The current for- 
mulation is ideal for this purpose, as will be illustrated 
later. 

Although the strategy is identical, there are subtle dif- 
ferences in the introduction of the formalism for one- 
dimensional systems, where we transform from the scalar 
wave function to a scalar variable, and higher dimensions, 
where the transformation is from the scalar wave func- 
tion to a vector variable. For this reason, we first intro- 
duce the formalism for one-dimensional systems and then 



afterwards consider the more general theory in higher di- 
mensions. 



A. One dimensional formalism 

In one dimension, we introduce the following real val- 
ued field related to the logarithmic derivative of the wave 
function, 



(3) 



where we use the notation that <p = Qipfyi unless the 
argument is specifically stated. 

The periodic boundary conditions on the wave function 
implies that <p is also periodic, but cannot be a constant 
function. Using <fi in the Schrodinger equation, we obtain 
the first order Ricatti equation, 



V{x) =E + <f>' + (j) 2 



(4) 



where we work in units of . Note that in these units, 

2m ' 

V and E have the dimensions (length) _2 , while <f> bas the 
dimensions (length) . 

Let us now consider (||). In principle, if we knew the 
functional dependence of the observable O on V, we could 
compute the desired quantities directly from (^|). How- 
ever, except in extremely trivial cases, we do not know 
this dependence and, in particular, we do not know the 
functional dependence of eigenvalues and eigenfunctions 
on V, which are the averages we would like to compute. 
On the other hand, the functional dependence of these, 
and many other observables, o n th e variable <p is fairly 
easy to determine (see section |l D| ). It therefore seems 
like a good strategy to change integration variables in 
(||) from V to <f> using the simple relation (Q). Doing 
so will shift the complexity of (g) from the observables 
to the action (probability distribution) for cj>, which will 
generally be highly non-linear, even though P[V] may be 
simple, i.e. Gaussian. The latter problem is, however, 
more amenable to treatment through the arsenal of per- 
turbative and non-perturbative field theoretic techniques 
than the original problem as stated in (0). 

To facilitate the change of variables (§) in (g), we in- 
troduce the identity^! 



1 = N~ 



= N- 



[d<j)}dES[ / dx(j)(x) ]6[V-E-(j)' - (j) 2 \\J\ 



1 l[d<t>]dE S[V - 



E- 



2 }\J\, 



(5) 



where the functional integral is over all possible non- 
constant periodic configurations of </>, N is the total num- 
ber of states (dimension of the Hilbert space), and the 
Jacobian | J| = | det(-^ + 2(j>)\. Note that the role of E 
under our change of variables is to replace the integration 
over the constant mode of V, which cannot be done with 
4>. It is easy to see, using the conditions imposed on <j>, 
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that the operator — ^ + 2(j> can be transformed to — -4- 
through a similarity transformation, thus not affecting 
the determinant. Therefore, \ J\ is simply a multiplicative 
constant, which can be combined with the normalisation 
of the functional integral. 

Inserting the identity (||) into the averages of (||), and 
then completing the integration over the disorder, allows 
us to obtain a field theory, formulated in terms of the 
variable 0, for the disordered average of observables 



(O) = Z- 1 J [dV]0[V]P[V] 

= Z- 1 J [dV] J-[d<p]dE 8[V-E-<j)' - <j> 2 ] d[V]P[V] 
Z- 1 l\d<p]dEd[E,(j)}P[E + (j)' + (f) 2 }, (6a) 



where 



}dE P[E + cj)' + cj) 2 ] = NZ/\J\. 



(6b) 



It should be noted that although we considered only a 
random potential when constructing this field theory, it is 
also possible to obtain the field theory when there is both 
a random potential V(x) and a deterministic potential 
W{x). In this case the result is similar to (pa]), except 
the energy is now shifted by W(x), 

(O) = Z- 1 l[d<p]dE d[E,<j)}P[E - W + </>' + <j> 2 }. (7) 



B. Higher dimensions 

In higher dimensions, we introduce a real valued vector 
field related to the gradient of the logarithm of the wave 
function, 
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(8) 



Since the wave functions are assumed to be of class C 2 , 
we note by direct computation that V x A — 0. Also, 
periodicity of tp demands that A is periodic and does not 
contain a constant mode. 

Using A in the higher dimensional Schrodinger equa- 
tion, we obtain 



V{x) = E + V • A + A- A. 



(9) 



As in the one-dimensional case, we can introduce an 
identity to implement a change of variables, based on (|^), 
between V{x) and the field A, constrained as described 
above: 



1 = N~ 



- 1 J-[dA]dE5[V - E - V ■ A- A- A]S[V x A]\J\, 

(10) 



with f- having the same meaning as in (S) . 

Using the identity (|l0|) in (^) and integrating over the 
disorder gives the corresponding field theory for the dis- 
ordered average of observables in higher dimensions, 

(6) = Z- 1 J[dA]dEO[E,A]\J\5[V x A] 

xP[E + V • A + A- A], (11a) 

with 

Z= J[dA}dEP[E + V -A + A-A]S[V X A]\J\. (lib) 

We can of course solve the constraint V x A = by 
setting A — V</> so that the resulting theory can be ex- 
pressed as a scaler theory. However, the form ([ll]) is 
particularly useful for the inclusion of a magnetic field 
B, since one only has to treat A as a complex field with 
the constraint V xA = replaced by V x A = ieB/c. We 
do not work out the details of this effective field theory in 
this paper, but rather focus on the one-dimensional case 
to illustrate the basic ideas. 



C. Translational Invariance 

Central to our analysis will be the translational in- 
variance of the effective action, which stems from the 
assumed translational invariance of the probability dis- 
tribution, P[V], and implies that translational invariance 
is restored after averaging over the disorder. The trans- 
lational invariance of the action leads to the appearance 
of an implied integration over a collective coordinate in 
the functional integral, corresponding to integration over 
the moduli space associated with the translational sym- 
metry. It is appropriate to make this integration over 
the collective coordinate explicit to ensure that is cor- 
rectly taken into account. To do this, we use a method 
inspired by the Faddeev-PopovEll quantization method of 
gauge theories. For simplicity we consider here the one- 
dimensional case, and indicate below how to extend to 
higher dimensions. 

We introduce the identity 



L/2 

1 = c / dx Q F'[<p Xn }S(F[(j) X0 } 

L/2 



(12) 



where <j) x ° = 4>(x + xq), F is an arbitrary functional of 
<p which is not translationally invariant, and F'[(f> x °] de- 
notes the derivative with respect to xq- Unless explicitly 
stated, the integration is over the interval [—•§-, §]• The 
proportionality constant c can in some cases be diver- 
gent, due to the existence of the Gribov ambiguitjO for 
certain choices of F. Under these conditions, one needs 
to be careful to extract the spurious divergent term af- 
ter completing the functional integration, as this term 
cancels with a similar term in the normalization. 
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As is usually done in gauge theories, it is more con- 
venient to implement the identity (|l2| ) after integrating 
both sides with c _1 J dvf(v), where / is an arbitrary 
function, so that 

c- 1 Jdvf{v) = J dx F'[ct> x °}f(F[n, (13) 

where the left hand side is independent of <fi. Multiplying 
both the numerator and denominator of (^aj) by the left 
hand side of (|l3| ) gives 



(6) = Zp l J-[d<P}dEdx F'[4> x °]f(F[<j> x °]) 

xd[E,4>]P[E + + <j) 2 ], (14a) 
where the partition function is given by 



Z F = j-[d4>]dEdx F , [4F°]f(F[,lf°°]) 

xP[E + cj)' + (j) 2 }. (14b) 

Note that the functional F[<j>] acts in a similar fashion 
to gauge fixing terms in conventional gauge field theories. 
Since the choice of F[4>] is arbitrary, we can choose F[<j>] 
so that our calculations can be simplified. As we shall 
see later, this choice will depend on what observables we 
wish to average. Also, any proportionality constants that 
appear due to the use of the Faddeev-Popov quantization 
method cancels out since we use the identity in both the 
numerator and the denominator (although, as mentioned 
earlier, extra care is needed if Gribov copying occurs). 

To extend to higher dimensions, we introduce a vector- 
valued functional F^ ] and the associated Faddeev- 
Popov determinant A[A] = det dl g i ^ . The results ob- 
tained above can then be generalized to higher dimen- 
sions through the replacement F^ ] — ► F^A 310 ] and 
F'[<f>] — > A[A], the latter denoting the corresponding Ja- 
cobian. 



D. Disordered averaged observables. 

Instead of writing the observables as a functional of V, 
we would like to obtain them directly as functionals of 
the fields (j> or A and the energy, E. It is possible to do 
this for observables like the density of states, correlators 
of the wave function, and the conductivity. 



1. Density of States 
The density of states at energy E is defined as 



where E m are eigenvalues of the Schrodinger equation, 
and d denotes the dimension of a system of size L. How- 
ever, if we consider the identities, (g, [n]), we see that the 
functional integral can be considered as a sum over all 
possible solutions of the Ricatti equations (^, ^), i.e all 
possible eigenstates with all possible eigenvalues for the 
corresponding Schrodinger equation, and is thus the total 
number of states. However, when fixing the integral over 
the energy at E, the functional integral yields only the 
number of eigenstates at E, and is therefore proportional 
to the density of states, p{E) for a particular configura- 
tion of the disorder. After integrating over V, this yields 
the disordered averaged density of states. Thus, we ob- 
tain via inspection the formula for the average density of 
states, normalised by the total number of states, for one 
dimensional systems, 



(P(E)) 



±fmp[E- 



(16a) 



and for higher dimensional systems 
1 



(P(E)} = 



ZL d 



[dA]P[E + V • A + A ■ A]5[V x A]\J\. 

(16b) 



2. Correlators of the wave function 

To obtain the observable related to correlations of the 
wave function, ip(x)ip(y), we use the definition of the field 
d>, (0) to write the unnormalised wave function as 



tp(x) = exp 



dx 4>(x) 



(17) 



Up to a global phase factor, all information that can 
be obtained from the wave function can also be obtained 
from (|l7j), including information about the phase, which 
we need to consider carefully when computing the local- 
ization length so as to avoid obtaining incorrect results 
due to random phase cancellation. 

The phase of the wave function changes as the wave 
function changes sign. From (|^), we see that <f> must 
diverge at these points. We thus need a prescription to 
calculate the integral in the exponent of ( |l7j ) at these 
points, since the result must be finite. The prescription 
we use is to integrate over a contour from to a;, where 
the contour avoids the positions on the real axis where 
there are singularities in <p(x) by moving around them in 
the upper complex plane with a semi-circle of radius e. 

This contour integral can be written in terms of the 
principle value of the integral plus a phase which depends 
on the number of times a singularity occurs in the interval 
[0, x\. Using this prescription, we are able to separate the 
phase from the integral over <p. Thus we obtain for the 
normalised wave function (the notation <f is used for the 
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contour as described above) 



ip(x) = y/N[4>] exp f dx'(j){x') 
Uo 

pf dx'<j>{x') -i7rVRes[ 
Jo 



N exp 



(18) 



where Xj 8X6 the positions of the singularities and 



dx exp 



2P/ dx'(j>{x') 
o 



(19) 



To avoid the problems associated with the random 
phase cancellation when computing the localization 
length, we need to calculate correlators between, the ab- 
solute values of the normalised wave functionsta. Using 
@H) with Q, we obtain an equation for calculating the 
disordered average of the 2-point correlator of the wave 
function at fixed energy E, 

\Mv)\) = 

z 



f-[dtj>] dx a F'[<j) X0 ]f(F[dp X0 ])N[dp{x)]P[E + <j}' + <jy 2 } 

(20) 











x exp 


pf dx'(j){x') 


exp 


pf dx'<f>(x') 




. Jo 




. Jo 



It is easy to check that this correlator is translationally 
invariant and that it only depends on \x — y\. 

Since F[cf) Xo } and / are arbitrary, we make a 
choice which simplifies the calculation of the correla- 
tor by cancelling out the normalisation factor, N[(jj\. 



To do this, we choose / = 1, and 
— J dz exp (2P f*dx'(f)(x' + xq)) , so that 



F 



'[<f°] = 2<j)(x ) Jdzexp 2P^ dx'<j){x' + x ) 



2<p(xo) exp 



o 

2P/ dx'4>{x) 

x 



Nttix)]- 1 . (21) 



Using ( pl|) in (20) we have 

(hM*)| \Mv)\) = 



Zp 1 l[d4>)dx (j){x )P[E + <f>' + <f> 2 







ry 1 


x exp 


pf dx'<t>{x') 


+ P dx'cj>{x') 









(22) 



Note that the integrand in ( p2| ) can be written as a to- 
tal derivative to xq, and thus the integral over xq naively 
gives a result of zero. This is, however, an artifact of the 
choice of the gauge in the Faddeev-Popov identity, which 
is zero due to Gribov copying. 

To obtain the correct result for the averaged observ- 
able, it is necessary, after completing the functional in- 
tegral, to extract the terms that give zero using some 
form of regularization, and cancel them out with similar 
terms that occur in the normalization. What remains is 
the correct result for the disordered average of the ob- 
servable. 

Using the periodic boundary conditions of 0, we see 
that we can transform ( p2] ) into a form where the sym- 
metries of the system are more explicit, 



(IMaOl \MV)\) 



dx K[4> X0 ] exp 

dx K[(j) X0 ] exp 
I — 



where K[(f> x °] = <f>(x )P[E + (/)' + 2 ]. From (|23^ ) we 
s ee t hat the correlator is translationally invariant, while 
( |23b| ) shows that the correlator contains a reflection sym- 
metry around \x — y\ — Thus one needs to ensure that 
any approximations that are made respect these symme- 
tries. 

These considerations can be generalized to higher di- 
mensions, with the integral in the exponent of ( |l8| ) being 
replaced by the line integral §~ A-ds along any path con- 
necting and x. Due to the constraint V x A = this 
integral is path independent. 

We note that if the wave function changes its sign (A 
becomes singular) along a certain path connecting and 





Ax- 


-v\ 


p 




dx 












\x-y 


p 


r 






Xq 





dx 1 <f>{x') + P I dx'(j)(x' 

xq 



(23a) 
(23b) 



x, it must do so along any other path, which implies 
that the associated singularity in A must appear in all 
possible paths connecting and x. This in turn implies 
that the singularity in A occurs on a surface separating 
and x into disconnected regions. Any path connecting 
and x may therefore cross a singularity and a prescrip- 
tion to handle this singularity is required. We can do 
this in the same way as the one-dimensional case : if t 
parameterizes the path, we can avoid the singularity by 
a detour in the complex plane. In this way the absolute 
value and random phase of the wave function can again 
be separated, with the principle value of the line integral 
determining the absolute value. The normalisation of the 
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wave function can again be cancelled by an appropriate 
(non-unique) choice of F so that the correlation in higher 
dimensions is given by 

(|^(x)| \i, E {y)\) = J [dA]dx A[A s °]f(F[A s °]) 
xP[E + V ■ A + A- A]\J\5[V x A] 



where 



x exp 



P\ A{x')-dx' + P A(x')-dx" 

xo J xo 



3. Conductivity 



(24) 



The conductivity of a system ofi— non-interacting 
fermions is given by the Kubo formulaE3. For our pur- 
poses it is convenient to integrate the Kubo formula by 
parts and use the periodic boundary conditions on the 
wave functions, ,so that the real part of the conductivity 
is then given bya 



Re<7(w) 



dE^±<S>(E,u;), (25) 



4e^ 



, d^ a 

dx 



E 

xS(E - E a )8(E + hcu - Ep). (26) 



Here f(E) is the Fermi- function, and all energies are mea- 
sured in units of Without loss of generality, we can 
focus on the disorder average of the quantity $(£ , ,w), 
which is of course just the contribution to the conductiv- 
ity of a particle at energy E. 

We are able to obtain the form of the observable for 
the disorder average of <& by using a technique similar to 
the one used to obtain the expression for the density of 
states. The only difference is that there is a summation 
over two eigenvalues in ( |2q), and so we need to introduce 
two identities of the forrnflq), where the integration over 
cj> is replaced by integrations over ef> a and <frp which corre- 
spond to two solutions of the Schrddinger equation with 
the same random potential. We can write the average of 
$ as 



where 



z- 1 Ud<t> a ][d<t> }5[hLu - 4>' a -<f>l + ^ + <P}]P[E + <j>' a + 4> 2 a ] 



dx'<fi a (x') + f dx'<j)p(x') 



4e^ 

hi' 



xJV^JJV^] / dx <j> a {x) exp 



dEdEj- [#«] [typ] S[E - 4»' a -€ + 4>'p + tf>]P[E + 4>' a + <t>l\ 



Completing the integral over <j)p, we have 

4e 2 



($(E,w)) = ^j-Z- 1 J-id^dxdxPiE + ^+^N^Ni^M^Mx) 

x exp dx'[4> a {x') + cf>p(x')]j exp (j> dx'[4> a (x') + 4>p(x')\ j . 



(27a) 



(27b) 



(28) 



r 



where 4>p(x) is a functional of <j) a and is determined by 

(29) 

We can now use the Faddeev-Popov method, where we 



Twj 



introduce the same choice of gauge for the <j) a and (j)p 
fields as in the previous section, which allow us to cancel 
out the N[cf> a ] and N[cf>p] normalisation factors respec- 
tively, so that 



4 e 2 r 

($(E,u)) = j^ z f X j-[d(t) a }dxdxdxodx <j) a (x)<j) a (x)(l3 a (xo)(t)p(x )P[E + (f>' a + <t> 2 a ] 



x exp 



dx'4> a (x') + <x> dx'cf>p(x') + <x> dx'<f> a (x') + &> dx'4>p{x') 

x J xq J xo J xq 



(30a) 
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with 



Z F = JdEdE l[d4> a }dx dx P[E + 0' Q + <t>l]F'[<Pl°]F'$f}. 



(30b) 



In higher dimensions we use the same strategy as above to obtain 



4e 2 

= hL 2_ d Z p l I [dA a ] dxdydxhdyhA a (x )A a (yo)A[A a ]P[E + 'V ■ A a + A a ■ ^4 Q ]| J Q | 



xA a (x)A a (y)exp 



A a (x ) -dx + & Apix) -dx + <b A a (x) ■ dx + j> A^x 1 ) ■ dx 
' Jyo J So Jyo 



(31) 



r 



where Ap = Ap[A a ) is the solution of 

hoj - V ■ A a - A a ■ A a + V ■ Ap + A f3 ■ Ap = 0. (32) 



III. ONE DIMENSIONAL SYSTEMS WITH 
GAUSSIAN DISORDER 

In this section we consider one dimensional Gaussian 
disordered systems. We do so to illustrate how the 
formalism as described in the previous section can be 
applied, using standard approximation schetp^, to re- 
cover known-results for the density of statesBE 3 ! and the 
conductivity!! 

If we have Gaussian disorder, then (V(x)V(y)) = 
l~ 1 S(x — y) and the probability distribution P[V] is given 
by 



Weak disorder limit 



P[V] = exp (-1 



L/2 



dxV{x) 2 ) 



(33) 



-L/2 



where the dimension oil is (length) 3 . Normally, as in the 
previous section, we are interested in observables at fixed 
energy, E. In this section we concentrate on these and 
therefore set 6\E, 0] = 0[E, 4>]5(E - E) in what follows. 
Using ( |33| ) in (g), or its equivalent form (|lj), we obtain 
a ip 4 field theory for calculating the disorder averages of 
fixed energy observables: 



(O) = z- 



\0[E,cf>}exp(-S[<f > ,E]), 



where the action is given by 

■ L/2 



S[(/>,E} =1 f dx(E + (j)' + dj?f 

J -L/2 



(34a) 



(34b) 



Here Z and O denotes the normalization and observable 



in a generic gauge and we omit the subscript F of (14). 

It is not possible to calculate the functional integral in 
( |34| ) exactly, so we use perturbative approximations in 
order to calculate the disordered averages. If / is large, 
equivalent to a weak disorder system, we expand around 
a saddle point in (0). For small (, or a strong disor- 
der system, we use a Hubbard-Stratonovitch transforma- 
tion to obtain a functional integral which can be approx- 
imated well in this regime. 



For largeLZI I, we calculate (34) perturbatively to one 
loop order using a saddle point approximation. Following 
a procedure similar to the one used in ^inn- Justing, we 
find approximate saddle point solutionsEJE3 which satisfy 
the saddle point equation to leading order in L^ 1 , and 
thus become exact in the thermodynamic limit. These 
saddle point solutions for positive energies are 



WW 



tan 



WE > 



(35a) 



where m is the nearest even integer to yEL/n, and for 
negative energies 



±V|-E , |tanhVl^l(^- *o) 

Kta,vh^/\E\(x + x ), V£<0(35b) 

where we assume a periodic continuation of (f>f (x) outside 
[—■j, —], and where the relative seperation 2xq hetween 
the instanton and the anti-instanton pair is largec3. We 
shall see below that the constraint that <f> contains no 
constant mode implies 2xq = — so that this condition 
is automatically fulfilled. Under these conditions of well 
seperated instanton and aitti-instanton pairs the dilute 
gas approximation is validcj. 

We wish to emphasize that the two solutions obtained 
for E > and E < have different physical behaviour, 
which can be explained by noting that the potential 
changes from parabolic for positive energies to a double 
well potential for negative energies. 

For positive energies, the requirement that the solu- 
tion must have periodic boundary conditions, as well as 
the constraint that there cannot be a constant solution, 



implies that the energies are quantized, leading to ( 35a ). 
Also, the condition that m is an even integer is due to 
the periodic boundary conditions of the Schrodinger wave 
function. 

We note that the solutions (|35a| ) are topologically dif- 
ferent for different m, as each solution has a different 
number of singularities. Since the number of singularities 
are related to the number of nodes of the wave function, 
each saddle-point solution corresponds to wave functions 
with a different number of nodes. To be more precise, 



8 



when compactifying the real line to S 1 by identifying 
+00 and — 00, we note that is a mapping from S 1 to 
S , where the mappings are classified according to wind- 
ing numbers. It is simple to see that the solution c/>™ has 
winding number m. 

In order to obtain a reliable approximation to the func- 
tional integral, it is necessary to sum over all the topo- 
logically different sectors. Additionally, since the fluc- 
tuations, 77, are smoothly varying functions around the 
saddle point solutions, we see that all the information 
about the phase (which is determined by the number of 
singularities in <ft) is contained in the classical solutions, 
so that the localization length (which is extracted from 
the absolute value of the wave function) is purely deter- 



mined by the fluctuations, rj. 

There is, however, a complication in the saddle point 
approximation, since the approximation contains a zero 
mode, which is a manifestation of the translational in- 
variance of the action ( |34b| ). To circumvent the prob- 
lems associated with the zero mode, we make the 
translational sym metry . explicit in the form of a col- 
lective coordinater 7 rr 8 H 29 l which we introduce by using 
the Faddeev-Popov method discussed earlier. Insert- 
ing the identity @ into (Q with the choice F[(j> Xo ] = 
J dxcf)o(x)(j)(x + Xq) in order to project out the zero mode 

cj)Q = ^jj-, we obtain the positive energy saddle point 
appoximation 



(6) « Z-'J2 [dx l[d V ]e X p(-lL(AE) 2 )F'[^+v]KF[€ l + v]) 

rn 

r L / 2 

xO[E,^ l (x-x ),r](x- x )]exp(-l dxTjA' 1 ^, \fE > 

J-L/2 



(36a) 



r 



where AE = E — ( m L ^) 2 and the propagator is given by 



A" 1 = (-— + 2£ + 6C 2 )- (36b) 



becomes the constraint that the rj integration is over all 
non-constant modes. Also note that in the pure limit 
(I — > 00) the value of m, and thus the topological sector 
is fixed by the energy. 



Here we noted that the constraint 5 (J dx(f)™) is trivially Integrating over the zero mode and making the approx- 



satisfied so that the constraint S( J dxcj)) in (34) simply imation that F'[<fi™ + 77] w F'lcf)" 1 ], we obtain 
I 



(O) 



Z- 1 ^ J " dx o Udq]' exp(-lL(AE) 2 )F'[ 



xO[E,<f>™(x-x ),r](x-x )]exp(-l / dxrjA^rj), VE > 

-L/2 



L/2 



(37) 



r 



where the accent denotes that the zero mode is excluded 
when calculating the funtional integral. Also note that 
F'[<fi™] is a divergent constant (the zero- mode is not nor- 
malizable), that may depend on m. However, requiring 
that the disorder average gives the correct result in the 
pure limit, we find that -F'[</>™] is a constant independent 
of m, which can thus be incorporated into the normal- 
ization. 

In the negative energy region, the saddle point equa- 
tion has a double well potential. The constraint that cf> 
has no constant mode does not allow us to obtain con- 
stant saddle point solutions situated at the minima of 
the potential, thiis | _t]jic ]| onljf r C(^lijer possible solutions are 
instanton solutionsc 



v,-athjer possible solutions are 
where tunnelling occurs 



from one minima to the other. Since we must also satisfy 
periodicity, there must also be tunnelling back to the 
original minima. This to and fro tunnelling can occur 
multiple timesE3, corresponding to topologically differ- 
ent sectors over which one has to sum, but since there 
is an exponential decay associated with each tunnelling 
process we consider only solutions, ( |35bD , where the tun- 
nelling occurs once. 

Once again the saddle point approximation contains 
a zero mode, which needs to be in tegra ted out. Addi- 
tionally, the saddle point, solution ( |35b| ) allows another 
quasi-symmetry to existEH. This quasi-symmetry is due 
to the fact that for large system sizes, local translations 
are possible that only changes the action by terms of or- 
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der e xp(— cL). To be specific, for large separations in 
( |35h| ), a translation in xq has an exponentially small ef- 
fect on the action, so that to leading order in i- it is 
a symmetry of the action, which in the L — > oo limit 
becomes an exact symmetry. Associated with this ap- 
proximate symmetry there is again an approximate zero- 
mode. 

To circumvent the problems associated with the zero 
modes, we make the translationa 
the form of collective coordinate! 
duce by using the Faddeev-Popov method discussed ear- 
lier. However, since there are there are two collective 
coordinates which we wish to introduce simultaneously, 
we need to modify the identity in ( |l2| ) so thatEZl 

,L/2 

l = c dxodxoAlcjylSiFl^^SiFl^ ^ }) (38a) 

J-L/2 



etries explicit in 
9 ' which we intro- 



where 



OF dF 

dF dF 

Oxq 3xq 



(38b) 



In order to project out the zero modes, we choose 

F[(j) xo ^ ] = [ dx(t> (x,x )<j)(x + xo) (39a) 
F\,.- '1 = I dxMx,x )cj)(x + x ) (39b) 



where the zero mode 4>o and quasi zero-mode <po are given 
by 



<j> (x,x ) 



Ox 



4> {x,x ) = ^. (39c) 
dx 



Using ( |38| ) and ( p9|) in (|34|), changing the variables </>(x + 
xq) — ► 4>(x) in the functional integral and then writing 
4> = (f>f + rj, we have 



(O) = 



Z- 1 Jdxodxo l[drf\d[E, 



4>f(x - x ) + rj(x - x )} exp (-lS[(f>f + r],E}) 



x A[0± + V ]S(F[^ + V })S(Flcj)t + 77]) 

I 



dx 



dx Q 



6{x - L/A). 



(40) 



Here we have handled the constraint 5(J dx((f>^ + rj)) 
by restricting the rj integration to be over non- 
constant modes, leaving the constraint S(J dxcf>f), 
which can explicitely be written as 5(J dx<p^r) = 

S(xq — L/A)/\ fdx-^-\. Note that the constraint on 
xq forces the instanton and anti-instanton pair of ( |35b| ) 
to be separated by -| , which for large L allows us to use 
the dilute gas approximations, 



(d)^Z- 1 \E\cxp(-^l\E\ 3 / 2 ) 



x J-[dri]'dx 0[E, (f)f(x - x ),r](x - x )] 

l-L/2 

xexp(-M dxrjA^r]) VE < 0, (42a) 

J-L/2 



where 



<t>f{x) w ±^/\E\taah.y/\E\{x- j)@{x) 

L, 



TVl £; l t anh v /|£;|(j;+ -)B(-x) VE <(61) 



After integrating c 
the dilute gas results 



2r xq in equation ( |40| ) 
1 one finds S[(j>f,E] = 



and using 

f\E\*/*, 



oc 



\E\ 3 / 2 and | Jdi 



dx 



OC 



the saddle-point approximation for the negative energies 



|.E| 1//2 , which gives 



A ± 1 = (-— + 2E + 6^t), (42b) 

and the notation [drj]' denotes that the zero modes are 
excluded in the functional integral. 



To calculate the functional integrals in (36) and fl42j), 
we need to be able to calculate the propagator, A, and 
the determinant of it's inverse, det(A _1 ). This involves 
solving the eigenvalue equation 



(43) 



where <p c is given by (35a) in the positive energy region, 
or by ([IT]) in the negative energy region. Although it is 
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possible to solve for the eigenvalues and eigenfunctions 
of ([l3]) exactly for the different energy regions, it is not 
possible to obtain closed expressions which is necessary 
when calculating the propagator. We thus need a con- 
sistent approximation for calculating the eigenvalues and 
eigenfunctions of (|4^) for both positive and negative en- 
ergies. (These approximations, as well as comparisons 
to the exact results, are discussed in detail in the Ap- 
pendix.) 

For positive energies, we note that the dominant prop- 
erty of the classical solution ( |35a| ) is that it contains 
singularities which appear periodically with period — . 
Thus, the eigenfunctions of p^ ) must be zero at the 
points where these singularities occur. Additionally, the 
periodicity of the classical solution leads to Bloch charac- 
teristics of the eigenfunctions, so that the eigenfunctions 
are also periodic over the interval which implies that 
we only need to solve (^) over the interval [— ^ , 5^] . 

We wish to make an approximation for the propagator 
that captures the essential characteristics of the original 
eigenfunctions. As a first approximation, we treat the 
tan 2 potential in ([l3]) to lowest order in perturbation 
theory, where we capture the singularities of the potential 
by imposing vanishing boundary conditions at We 
thus need to solve the eigenvalue equation 



A„*„, 



(44) 



with vanishing boundary conditions at ± We can now 
calculate the approximate propagator and determinant in 
the positive energy region to give 



A m (x,y) 



dx 2 



2E 



2 ~ /cos([2n- l]2^a;)cos([2n- l]^y) 



L 



([2n- W- 
sin(2n^a;) sin(2n^%) 



2E 



(2«T) 2 



and 



|det(A" 1 ) 



2E 



■ sinh 



*l J 



(45a) 



(45b) 



In the negative energy region one finds from the exact 
solutions two zero modes (which we must eliminate) and 
a doubly degenerate bound state at A = 3|1£|. The rest of 
the spectrum consists of four fold degenerate scattering 
states starting at A = &\E\. From the exact solution it 
turns out (See the Appendix) that the eigenvalues and 
eigenfunctions of the scattering states can be well ap- 
proximated by that of a free particle. This amounts to 
writing tanh 2 (x ± xq) = 1 — sech 2 (a; ± xq) in the clas- 
sical solution appearing in ( }42b| ) and then treating the 
sech 2 (x ± xo) terms to lowest order in perturbation the- 
ory. Thus the approximate eigenvalue problem that has 



to be solved to obtain the spectrum and eigenfunctions 
of the scattering states is 



-—+4\E\)V n = A„*„, 



(46) 



where ^ n are periodic on a system of length L. Tak- 
ing into account the bound states and scattering states, 
approximated as a free particle spectrum, we obtain the 
determinant in the negative energy region as 



|det(Aj 



oc sinh 



(47a) 



The contribution of the bound states in the propagator 
can be neglected as it involves the product of two very 
well localized eigenfunctions evaluated a points which are 
well seperated. Taking into account only the scattering 
states, again approximated as free particle states, we find 
for the propagator 



A±(x,y) 



E 



cos(2p(x-y)) 

(2^)2+41^1 



(47b) 



Using these approximate results for the propagator and 
determinant for positive energies (^5|) and negative ener- 
gies (|47j) with the respective formulations for disordered 
averages, (|3^) and (j4|), now allows us to calculate the 
average value of observables to leading order in L. 

As mentioned earlier, this approximation is valid when 
I is large in a sense determined by the other two scales, 
namely E and L. To find the precise criterion one has 
to evaluate the higher order loop corrections to ( |37| ) and 
([42|), e.g., in the case of the density of states, one has to 
evaluate higher order vacuum diagrams. Doing this, one 
finds that these contributions can be neglected under the 
condition that ^ > 1. 

There are two limits under which this condition can 
be fulfilled. Firstly, for a fixed energy E, we find that 
I ^> -§, which is large for large system sizes. The approx- 
imation thus holds in the weak disorder limit. Secondly, 
for a fixed amount of disorder, we have that E ^> h. We 
thus find that the approximation also holds in the high 
energy limit. Note, however, that in the thermodynamic 
limit (L — > 00), the condition cannot be satisfied, unless 
either I — * 00 or E — > 00, implying that using this ap- 
proximation for disordered systems only holds for finite 
system sizes. 



1. Density of states 

We now apply the saddle-point approximation dis- 
cussed above to the disorder averaged density of states 
given by ( jlq). For positive energies, this amounts to 
computing ( B7p with the observable = 1. Upon inte- 
grating over the functional integral to obtain the factor 
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det(A^ 1 )^ 1 / 2 , we have 

( P (E)) = N + J2^p(-lL[E-(^n 



2 EL . / V2EL\ 



be written as 



1 1/2 



-cosech 
\m\ \ \m 



(,48a) 



where N + is an unknown normalization factor, indepen- 
dent of the energy, that needs to be fixed in some manner. 

For negative energies, the saddle-point approximation 
of (|ilf ) leads to (^) with = 1. After integrating over 
the resulting functional integral, we obtain 



(p(E)) = N_\E\exp 



--l\E\^ 



cosech 2 (y/\E\L^ , 

(48b) 



where iV_ is another energy independent normaliza- 
tion factor. The latter result agrees with the result of 
HalpcrinEI obtained by different means, although his re- 
sult does not include the higher order corrections. It 
should be noted thatin order to compare with the re- 
sults of Lifshits et alc3 : which are given for a system of 
infinite size, one must take the L — > oo limit in the results 
above such that the ratio El/L is fixed. 



2. 2-point correlators 

If we consider the disorder average of the 2-point cor- 
relator given by , we note the exponential terms can 



exp 
exp 



P dx'4>(x')+P dx'<j){x') 



P dx'<f> c (x')+P dx'4> c (x') 



dx'f](x') + / dx'rj(x') 



x exp 



where we have relabled the the xq integration in (|22j ) 
to z. The linear terms of r; in the exponential can be 
written as an integral over the interval [-rS -§■] by using 
a combination of step functions, 



■v rL/2 
dx'r](x')+ I dx'r](x') — / dx'S(x, y, z\x')r](x'), 

J-L/2 

(49a) 



where S(x, y, z\x') can be considered as a source term for 
the rj fields and is given by 

S(x, y, z\x') = 9(x - x') + 6{y - x') - 26{z - x'). 

(49b) 

Furthermore, a periodic continuation is understood out- 
side [-<§,£]. 

Ap ply ing the positive energy saddle point approxima- 
tion ( p7\) to (|22|), and making the change of variables 
z — > z + xq , we have 



L /2 

(\^ E (x)\ \^ E {y)\) = Z- 1 ^ J dx dz^^[dv]'exp(-lL(AE) 2 )exp(-l J ^ dx^rj) 





r r° 










x exp 


2P dx'(j)(x') 


exp 


pf dx'[^(x' -x ) 


f 77] +p[ V dx'[cj>™(x' -x ) 


FT?] 








. J '-co 


Jx 





(50) 



Note that as a result of the gauge that we chose to can- 
cel the normalization of the wave functions giving rise to 
(p2|), equation (^0|) contains a total derivative with re- 
spect to z, which na'fvly gives a result of zero when com- 
pleting the integration over z. However, since we have 
a ratio of total derivatives, we should obtain a non-zero 
result if we use a consistent regularization method. With 
this in mind, we integrate over z and cancel the result 



with a similar term in the denominator of (J5CJ). We can 
now integrate over the fluctuations to obtain 

(We(x)\\iI>e(v)\) „ C(x,y) 

(\M0)\ |^(0)|) ~ C(0) 1 ) 

where 
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C(x,y) = ^ dx cos(J^-{x - x )j cos(J^-(y - x )j exp (-IL(AE) 2 ) 

m 

x(det(A- 1 ))" 1/2 exp 
After using (Bq) for the propagator and determinant, we have 



— jdx'dx" S(x, y, x \x')A m (x', x")S(x, y, x \x") 



where 



C(x,y) = ^2 J dx o cos(^-(x - zq)) cos(^-(y - x )^ exp(~lL(AE) 



'2EL 



cosech ( X * 



"I 1/2 



V h ;j cxp U F(wo) )' 



2 1 T ( i >rrnr. \ f mit 

F(x,y,x ) = j^}^J^ 2 — - [cos(J2n- 1) — (a-x )J +cos(J2ri- 1) — (y 



x ))-2 



2 1 r / mir \ ( run ' 2 

-— Yj— |sm^2n— — (x - x )J + sm[2n—(y - x Q 



(52) 



(53a) 



(53b) 



r 



and 



D n = 



L J 



2E 



/ TITYLTT \ ^ 



In the negative energy region, the saddle point approx- 
imation (E2) is obtained using the dilute gas approxima- 



tion. Ho weve r, we find that the approximate saddle point 
solution ( |35h| ) in a dilute gas approximation, (^l|) , breaks 
the symmetries of the system, thus we first need to write 
the exponential terms in (J22) in a form where the sym- 
metries are explicit. We do this by using (p3|) so that 







ry 1 


2 exp 


p[ dx'(f>{x') 


+ P dx'(f>{x') 




. Jx 


J x 



exp 



P 



x~y\ rO 




r\ x -v\ 


dx'<p(x') + P / dx <p(x) 


+ exp 


P 1 dx' 4>{x 


Jxq 




Jxq 



(55) 



Using this form for the exponential insertions in ( p2| ) , along with saddle point approximation for negative energies 
(E3), and then integrating over the fluctuations, gives 



C(x,y) = jdxoex-p 
I dxo exp 







\x-v\ 


p 


r 


( 




X 






r\ x - 


-y\ 


p 




dx' 




x 





>A±f„> 



Jx±tJ 



1 

4l J 



1 

4T 



where 



4 1 
F(x,x ) = 7^77 



n—l 



COS 



/2mr 



x ) — 2 cos 



2nn 



(x — .Tq) I — 2 cos 



2nn 



-x 



(56a) 



(56b) 



and 



Note that in obtaining (d9), we once again relabled the 



D n = 



2mr , 

— ) + 4|£| 



. ,57, 
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integration variable in (£2|) from xq to z, translated 
z —>■ z + xq, and then, as above, cancelled the z depen- 
dent terms giving rise to a total derivative with respect 
to z. Note from the saddle point solution ([ll]) that to 
leading order this correlation function decays or grows, 
like exp(±y/Ex), confirming the result of Lifshits et alc3 
that the localization length is proportional to 1/ \fE for 
large negative energies. 



and a fluctuating term, i.e. <f>p = <f>% + \- Using this 
ansatz, as well as the expansion 4> a = 4> c a + rj in (pjJ), 
and neglecting the coupling terms between the classical 
and fluctuating terms, we find that x = V an d <j>% must 
satisfy the saddle point equation with the energy shifted 
by huj, 

hu - (C) 2 - FJ + (<f>%) 2 + 4> c b = o, (58) 



3. Conductivity 

To calculate the disorder average of the conductivity 
given by ( J30| ) in a saddle point approximation, we first 
need to be able to solve for </>^ = 4>p\(j> a \ using (|29|). To do 
this, we make the ansatz that <f>p consists of a classical 



with the positive energy solution given by 

^(x) = -^tan(^ x). (59) 

To satisfy the periodic boundary conditions, we see that 
q must be the nearest integer to (m 2 + TiujL 2 jn 2 ) 1 ! 2 . 



I 

Applying the saddle point approximation in the positive energy region ( |37| ) (where we relabel xq to z) to (^), and 
making the change of variables x — > x + z, x — > x + z, xq — > xq + z and xq — > xq + z allows us to integrate over z 
(since the integrand is now independent of z) to obtain 



E m expHL(A£) 2 )<lH 



($(E,uj )) ^ m exp(-;L(AS) 2 )l>(a;o)' 



(60a) 



with 



$(w) 



d 

[dr]]' dxdxdxodxo — — exp 
dx 



dx'(^(x')+r,(x')) 



dx 



exp 



dx'(^(x')+r,(x')) 



x [<f%{x) + n(x)] exp 
x [4>a (x)+v{x)} exp 



dx'<j>%{x') + f dx'<f>™(x') + f dx'^(x') + f dx'^{x') 
ii Jo Jo Jo 

L2 dx' [l v &)A£r,(x>) - S(x')rj(x')} 

-L/2 

I 



(60b) 



where S(x') = S(x,x,Q\x') and the classical solutions in 
the positive energy region are now denoted by the su- 
perscript m. Note that the total derivatives that appear 
are due to the original Faddeev-Popov method used to 
cancel out the normalizations of the wave functions. As 



before, we can cancel them with similar terms in the de- 



nominator of (60a). 

Introducing source terms for the rj(x) and r](x) terms, 
integrating over the fluctuations and then integrating by 
parts, we have 



where 



F(x,x) 



dxdx exp [ jF{ x j x ) 



<2EL , / V2EL 
-. — |— cosech — : — — 
\m\ \ \m\ 



qn /qn \ (mn \ 

T sm {T x ) cos {— x ) 

qn / qn _\ (mm 

T SU \T X ) C0S {— 



mn /"ran \ ( qn \ ~ 

— sm(—x) co S [ T x)_ 

mn fmn_\ / qn _\~ 

— sm(—x) co S [ T x)_ 



2 
L 



El r . / mn \ . / 
sin \2n x)+sm[2n- 
, D 2n V \ L J \ 



n=l 

oo 

E 



L 



D 



2n-l 



I In 1 \ m7F \ (in 1 \ TO7r _\ 

cos^(2n — l)— — xj + cos(^(2n — l)—jr-xj 



(6la) 



(61b) 
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with D n given in (p4). 



with the observable given by 



B. Strong disorder limit 

In this limit, when I is small, we use a Hubbard- 
Stratonovitch transformation on (p4|), so that 

1 



(O) =Z- X /[dA]0[E,A]exp(-- \dx\A 2 + 2iElA]] 



(62a) 



The A dependent observable is given by 
0[E, A] = f\d0}O[E,0} 



x exp(— / dx 



I / + 2i(j) 2 A 



) (62b) 



and 



[dA]dE[det(-l 



dx 2 



x ex P(-y Jdx[A 2 



-21A)]- 1 / 2 
•2LE/A]). 



(62c) 



O(E,A )=j-[dcf,}d[E, 0] 



exp 



dx{l(j)' +2i0 2 A o ) 



(64b) 

Assuming that 0[-E, 0] is real, we note that 
O*(E,A ) — Q(E, — A ), which allows us to integrate 
over positive Ao in (64a) if we take the real part of the 
integrand, thus 



(O) = Z 



1 ' dA ReO(E,A ) 
o 



x exp 



.^-2i£A L 



(65) 



Here we considered only the lowest order approxima- 
tion where we totally neglected the contribution from 
the 4> 2 r\ term, but it is easy to extend d&jj ) to include 
the contribution of the quadratic terms in rj arising from 
the 4> 2 r) coupling, upon which the rj integral can still be 
done to yield a determinant, which will give higher order 
corrections to (p4|). 



We now approximate (^2|) by expanding (as described 
in Zinn-JustinEj) up to first order in the loop correc- 
tions. We do this by first splitting the integral over A 
into an integral over the constant mode, Ao, and non- 
constant modes rj. This is achieved by inserting the iden- 
tity J dAgdri5\A + A p\5( f dxrf) into the numerator and 
denominator of (|62a| ). Integrating over A, we have 



(O) 



dA Q l[drj\0[E, A + rf\ exp 



1 



dx rf 



x exp 



AjL 
I 



2iEA L 



(63) 



where the observable 0[E, Aq + rf\ is given by fl62b| ) with 
A = A + 77. 

We wish to approximate ( |63| ) by neglecting the 4> 2 r) 
coupling term that appears in 0[E, Aq + 77]. By exam- 
ining the effective action for Ao, expanded around the 
equilibrium value for Ao, we find that one can neglect 
the coupling term if ^ 1. As in the saddle point ap- 
proximation, this condition can be fulfilled in two limits, 
namely the strong disorder limit (I <C -g), or the low en- 
ergy limit (E j). Also, in the thermodynamic limit, 
this condition always holds for fixed energy and disorder. 

If we neglect the 4> 2r q coupling term, we can integrate 
out the rj integal, so that the disordered average is 



(O) = Z- 1 dA O(E,A )exp 



A 2 L 



2iEA L 



(64a) 



1. Density of States 



For the average density of states, the observable (p4q) 



is 



O(E,A ) 



det 



d 2 2iA„ 



oc 



dx 2 I 

(l+iyjoi 



2y/i 



cosech 



1-1/2 



(l + iysoi 



2 vi 



(66) 



which we can use in (pq) to obtain 



(p(E)) = / dA Re 
Jo 



[1 + i) cosech 



/ (l + i)VA^L \ 
I 2Vl J 



A 2 L 

x v/A exp ( 2i£A L 



2. 2-point Correlations 



(67) 



The observable used to calculate the correlator in the 
dual region can be obtained from (|2^) . Explicitly w ritin g 
the total derivative and using the approximation in ( |64lj ) , 
we integrate over the <fi field to obtain 
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O(E,A ) 



d (1 + i)^ /(l + i)V^i 

axo — - cosech ' 



x exp 



dx 
4/ 



2Vl \ 2\fl 

dx'dx" S(x, y, xq\x')A(x' , x")S(x, y, x$\x") 



(68) 



where A(ai,z") = (--£? + life.)- 1 . 



Using this in (p5|), extracting the terms that are xo 
independent and then cancelling the xq integral with a 
similar term in the normalization, we have 



\^ E (y)\) = C{x,y) 
(|^(0)||^(0)|) C(0,0) 



(69a) 



where 
C(x,y) 
x exp 



° d A Reii±^cosech^ 1 + i ^ L 



A 2 L 



2Vl 
2iEA L 



2^fl 



exp 



j t F{\x-y\ 



(69b) 



with 



F(x) 



cos(2^x) 



\(2HJL\ 2 _L 2iA 1 ( 2mr 
n=l [\ L ) + I \ \ L 



(69c) 



C. Microscopic realisation of the model 

It is useful to have a microscopic realisation of the 
model presented in the previous sections that relates the 
parameters to microscopic quantities. For this purpose, 
we consider a one dimensional model of N Dirac-delta 
scatterers placed randomly on a ring, with the potential 
given by 



V{x) = a^^8(x - Xi) - 

i=l 



(70) 



Here a is a dimensionful constant (units of (length) -1 ) 
that determines the strength of the scatters, and the sub- 
tracting term is chosen so that (V) = 0. Applying our 
formulism, the disordered average of some observable O 
is now given by 

.1 Jf , „ 

(6) = z- 1 j 2 L l[^ J-[d4>}dEd(E, 0) 



xS 



E 



2 ?:=i 
i2 



N 



a S(x — x , 



aN 
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(71) 



where 



,L- N 



z = j^Yi^-fmdE 



xS 



N 



E 



+ (f> 2 - S{x - Xi) + 



aN 



i=l 



(72) 



Introducing a Fourier representation for the functional 
Dirac-delta, we can write (0) as 



(O) = Z~ L j-[d4][dk)dEO{E^) exp 
x exp i / dxA { E 



AT log 
aN 



dx 



— iaA 



(73) 



with a similar expression for Z . Assuming that the scat- 
ters are weak (a small) and the system size L is large, we 
expand to lowest order in a and 1 / L : 



(6) = Z- 1 j-[d(f>}[dA}dEd(E,(j)) exp 



i / dxA E- 



Na 2 
~2L 



dxA* 



x exp 



(74) 



Performing the A integration, we have: 
(6) = Z- 1 i[d<t>}dEd{E,4>) 



x exp 



2Na 2 



dx (E + cj)' + (j) 2 ) 2 



with 



Z = j-[d(j)]dE exp 



(75) 



(76) 



We can therefore identify the disorder parameter, I, of 
our Gaussian model with I = L/2Na 2 — //2a 2 , with I 
the mean free path length. 



IV. NUMERICAL RESULTS 

In this section we numerically calculate the main re- 
sults that we obtained in the previous sections and gen- 
erate plots from these calculations in order to obtain a 
understanding of the results. 
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Density of Stales (L = 10) 




FIG. 1: Plot of the density of states (in arbitrary units) vs 
energy for L — 10. In (a), the result of the positive energy, 
weak disorder approximation for various values of I are shown. 
In (b), we plot the strong disorder approximation result, while 
(c) shows the result for the negative energy, weak disorder 
approximation. 



A. Density of States 

There are three different regions that we need to con- 
sider when calculating the disordered averaged density of 
states. Firstly, if E 3> y, then the weak disorder saddle 



point approximation in the positive energy region, (48a), 
holds. For E -C — j, the weak dis order saddle point ap- 
proximation for negative energies, (48b), is used. Finally, 
when — j <C E <C —, we use the strong disorder approx- 
imation, (|67|). Thus, we find that the weak disorder sad- 
dle point approximation describes the high energy tails 
of the average density of states, while the strong disor- 
der limit describes the low energy states. Note that for a 
fixed disorder, the only approximation that holds in the 
thermodynamic limit, L — > oo, is the strong disorder ap- 
proximation, as the regions described by the saddle point 
approximation tend to negative and positive infinity. 

We calculate the disordered averaged density of states 
for finite L, with various disorder values using the ap- 
proximations in their respective energy regions. These 
results are shown in Fig. [I] in arbritrary units. Fig- 



ure [j](a) shows the plot for the positive high energy tail 
for various disorder values. Note that at high energies, 
the density of states is peaked around the discrete values 
that one would get in the pure limit. Also, the width 
of the Gaussian distribution around these discrete values 
increase as the disorder in the system increases. Eventu- 
ally, when the disorder is large enough (or the energy low 
enough, as seen in the plot) the Gaussian distributions 
start to overlap, which leads to a change in the density of 
states from a almost pure behaviour to a strong disorder 
behaviour. Note, however, that the criterion that E ^> j 
no longer holds in this region, and that the strong disor- 
der approximation should be used instead. Figure |I](b) 
shows the density of states for the low energy states. As 
the disorder increases, the width of the strong disorder 
density of states increases. In Fig. 0(c), the result for the 
negative energy tail is shown. In this region the density 
of states falls off exponentially to zero. To describe the 
density of states over all energies, the arbitrary normal- 
ization factors appearing in the three different regions 
should be fixed by requiring a continuous matching at 
the transitional points E = ± j and by imposing some 
global normalization condition. Here we have only im- 
posed an arbitrary normalization within each region to 
exhibit the main features of the different regions. 

From the plots in Fig. [j], we see that there is a crossover 
from the almost pure system behaviour at large positive 
energies to an exponential decay at negative energies. Of 
particular note is the crossover region at small energies 
where the strong disorder approximation is valid lead- 
ing to a non-zero result for the density of states at zero 
energies. 

Also, we note that as the parameter I decreases, the 
width of the density of states in the negative energy re- 
gion increases due to the creation of additional bound 
states in the more disordered system, whereas for large 
I, that is a more pure system, there are less states at 
negative energies and the peak increases, leading to the 
E^ 1 / 2 singularity at zero energy for pure systems. 



B. 2 point correlators 

The correlation function in the weak disorder saddle 
point approximation is given by ( |5l| ) using (|5^) for pos- 
itive energies, and (|5^) for negative energies, while ( |69| ) 
gives the correlation function in the strong disorder ap- 
proximation. Without loss of generality, we can set y = 
and x = dL, allowing us to calculate the correlation func- 
tion in the appropriate energy regions where the various 
approximations hold. These results are shown in Fig. ^ 
and Fig. 3. 

Figure 2 shows the result of the correlation function at 
a fixed value of I and L for various energies. Figure ||(a) 
is the correlation function calculated in the weak disorder 
saddle point approximation with a large positive energy 
(E > £ ). In this region the result is dominated by the 
pure solution, giving rise to the oscillations. The effect of 
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Correlation Function (Length = 10) 
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FIG. 2: Plots of the correlation function (normalized to 
one at d = 0) vs distance d showing the disorder depen- 
dence, where (a) is the positive energy, weak disorder sad- 
dle point approximation (£7 — 1.5 x 10 3 ), (b) is the strong 
disorder saddle point approximation (£7 = —2.5) and (c) is 
the negative energy, weak disorder saddle point approxima- 
tion (E = —1.75 x 10 3 ). The disorder value in all three plots 
is / = 1.0 x 10 2 , and L = 10. 



the disorder is minimal and leads to a weak modulation 
of the envelope. Note that the solution is periodic, with 
a period shorter than dL, with only one period being 
plotted in Fig. |(a). 

As the energy is lowered even more, the region where 
the strong disorder approximation (— j <C E <C j) holds 
is reached as shown in Fig. ||(b). The correlation function 
decreases as the distance increases until d = 0.5L after 
which the correlation increases again. This is of course 
due to our ring topology. Since the bulk of the states 
occur in this region, we consider the behaviour of the 
strong disorder region in the next figure, Fig. ||. 

Figure ||(c) shows the correlator function for large neg- 
ative energies, where the weak disorder saddle approx- 
imation E <C —j once again holds. Here the results 
show that the correlation function decays exponentially, 
where the decay length is determined by the energy and 
is largely disorder independent. This exponential be- 
haviour is due to the formation of bound states in the 



Correlation Function (Length = 10) 
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FIG. 3: Plots of the correlation function (normalized to 
one at d = 0) vs distance d showing the energy and disorder 
dependence of the strong disorder approximation, where (a) 
the energy is fixed (£7 = 0), with various disorder values, and 
(b) the disorder is fixed (I = 2 x 10~ 2 ) and the energy is 
varied. The length of the system in both plots is L — 10. 



disordered potential. 

In Fig. ||, we consider the behaviour of the correlator 
in the strong disorder approximation, when the disorder 
and energy are varied. In Fig. |^(a), we keep the energy 
fixed and vary the disorder. As can be expected, there is 
a stronger decay when the disorder in the system is in- 
creased (smaller I). In Fig. ||(b), the disorder parameter 
is kept fixed, while the energy is varied. Once again, as 
is expected, the decay increases as the energy is lowered. 
Thus the strong disorder approximation gives a bridge 
from the almost pure behaviour at high positive energies 
to the strongly localized behaviour at strong negative en- 
ergies due to the formation of bound states. 



V. CONCLUSION 

In this paper we introduced a functional integration 
formalism for studying disordered averaged observables 
that provides a complementary viewpoint to the standard 
field theoretic techniques used at present. The formalism 
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is based on changing variables from the random potential 
describing the disordered system to a new set of random 
variables related to the logarithmic derivative of the the 
wave- function. This allows a more direct computation of 
certain disordered averages, such as the density of states 
or obsevables that explicitly depend on the wave func- 
tion. In particular we showed how to calculate the dis- 
order averages of the density of states (|l6|), the 2-point 
correlators of the wave-function ( 25 , |24|) , as well as the 
real part of the conductivity ( f30| , 31). 

As an illustration of how the formalism works, we con- 
sidered one dimensional Gaussian disordered systems. 
We were able to obtain results for the weak disorder and 
strong disorder limits for the density of states ( [48| , 67), 
and the 2-point correlators ( f33| , [S6| , |69|) . Unfortunately 
we were only able obtain results of the conductivity in the 
weak disorder limit ( |6l"l ) , as there is a complication in the 
perturbative expansion of the strong disorder limit when 
using the Hubbard-Stratonovitch transformation on (30), 
which we have as yet been unable to resolve. The formal- 
ism reproduced the results of HalperinEl for the density 
of states, and considering the the 2-point correlator we 
showed that in the thermodynamic limit all states in one 
dimension are localised. 

Future developments include the addition of a deter- 
ministic potential to the formalism, allowing magnetic in- 



teractions to be included. Also, the calculation in higher 
dimensions needs to be investigated further to see if signs 
of a metal insulator transition can be found. 
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APPENDIX A: PROPAGATORS AND 
DETERMINANTS 

We wish to calculate the propagator, A , an d determi- 
nant, det(A _1 ), where A -1 is given by (36b) or (421) 
This involves solving the differential equation (|43|), 



(Al) 



with periodic boundary conditions over the interval 

2 ' 2 J 



In the positive energy region <f> c is given by ( 35a ) , so 
that (|Al|) is 



d 2 ( ran \ 2 
^ +2 (-) ^ 



rrnr (vim 
— tan(— x 



• 71 



A n + 2(— J -2E 



(A2) 



The differential equation in (A2) can be solved exactly using the method of generalized ladder operatorsciH so 
that 



, t , /m7r\ 2 

A„ = (n* + 6n + 3) [—) + 2E Vn>0 



(A3a) 



and 



n-l 

k n 



cos I — —X 

d , , rrnr / mir 



i=0 



, +(i + 3)— tan 
ax L 



n+3 



V n > 0. 



(A3b) 



Unfortunately, although we have the exact solution, 
we do not have the eigenfunctions in a closed form. This 
makes the calculation of the propagator difficult. We 
thus wish to make an approximation to ( A2 ) that allows 
us to calculate the propagator. Since the tan 2 potential 



in (A2) contains m singularities in the interval [—-k, 4], 



the eigenfunctions must be zero where these singularities 
occur. The approximation that we make for (A2) must 
preserve this global property. The approximation that we 
make is to ignore the tan 2 term in (A2) and to change 



the boundary condition so that eigenvalues have zeros at 



the correct intervals. The eigenvalue equation that we 
must solve is thus 



d 2 

with the boundary condition that \I>(±^-) = 0. 

Excluding the constant mode as required by ' 
solution of (A4) is 



A n 



L ) 



2E V n > 



(A4) 
f), the 
(A5a) 
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with 



L /nmn 
= \ I — cos I - 



V L 



' L /nmn 



V odd n > 

V even n > 0. (A5b) 



As expected (A5a) and (A3a) arc in good agreement for 
large n. 

The propagator can now be calculated using the energy 
representation to give 



A m (x,y) = \~— + 2E 
2 ~ ^cos([2n - l]2Ht x) cos([2n - l]*f-y) 



sin(2n 



([2n-l]T) 2 H 
^x)sm(2n^y) 



2E 



2E 



(A6) 



while the determinant can be calculated with various 
techniques, i-e. the identity 1.143.1 found in Gradshtcyn 
and RyzhikEj, so that 



|det(A^)|=C- 



'2EL 



■ sinh 




(A7) 



where C is an energy independent constant, but not nec 
essarally independent of m. We de termine C by requir- 
ing that the density of states (48a) be the correct solu- 
tion in the pure limit. Obtaining the pure solution (up 



to a global normalization constant), requires that (A7) 



must be a constant independent of m. However, since 
E = (^jf-) in this limit, and C is independent of E, we 
hnd that C must also be independent of m. Thus we see 
that C is constant that is independent of E and m. 

In the negative energy region, C is given by (^) where 
we are us ing the dilute gas approximation. Thus the 
equation ( |Al| ) becomes 



dx 2 



•4LEI 



= A n *n, 



(A8a) 



where the quasi zero-mode 4>o is given by 



sech 2 
sech 2 



\E\(x-±) )B(,) 
E|,z + |))(-),-,S. 



(A8b) 



If we integrate over (A8a), and use the condition that 
</>o is orthogonal to we find that the constraint that 
the eigenfunction cannot contain a zero mode is satisfied 
for all eigenfunctions except the one corresponding to 
A n = 4\E\. The eigenfuction with this eigenvalue must 
be explicitely checked to see if the constraint is satisfied. 

We first calculate the eigenfunctions and eigenvalues 
of (A8) in the subinterval [— 4Jl]_or [0, |-] using the 
method of generalized operatorsEjO or via the solution 
of a hypergeometric equations, and find that the eigen- 
values consists of two discrete eigenvalues, Ao = and 
Ai = 3|i£|, and a continuum of eigenvalues A^ = 
(k 2 + 4)\E\. The corresponding eigenfunctions are 



cx sech* I y/\E\(xT t) ) ( ±a O> *ftx 



tanh{y\E\(xTk) 
cosh(^/\E\(xTT) 



Q(±x) 



(A9a) 



and 



3 tanh 



E\i.r?j)\ - I k 2 



^ kx e{±x) -3ifc tanhf y/\E\(xT j) ) e i ^ Flkx Q(±x). (A9b) 



We can now calculate the eigenfunctions for the full region [—-§,-§] by matching the eigenfunctions in (|A9] ) at x = 0, 



i.e. requiring that *+(0) = *~(0) and ^-(0) = ^-(0). The eigenfunctions are 
*o oc ±sech 2 ( \/\E\{x + -) ] Q(-x) + sech 2 ( y/\E~\( x ~ t) ) @ ( x ) 



*i oc ±- 



tanhf y|^[(a; 



cosh(VW|(a;+ §) 



9(-z) 



tanhf ^AE~\(x - h 



cosh 



jh(vPIOr-£)) 



e(x) 



(AlOa) 
(AlOb) 



* fc cx ±(2 - fc 2 )eV l B l fc (*+ 2 )0(-a!) ± 



30(-a;) 5 



\E\ dx 



+ (2- fc 2 )c i Vl £ l fc ^-7)e(x) + 



36(a;) d 



tanhf - ~) ) e 1 ^^-*) 



(AlOc) 
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where k ^ since the corresponding eigenfunction does 
not satisfy the constraint that there are no constant 
modes. Also, the periodicity requirements on the eigen- 
functions imply that k must satisfy the equation 



\E\kL 



2 - k 2 + 3ik 
2 - k 2 - 3ik 



(All) 



Note that there are degenerate solutions for each eigen- 
value, since the eigenfunctions can be constructed as a 
symmetric or an anti-symmetric solution. 

As in the positive energy region, we do not have the 
eigenvalues and subsequently also not the eigenfunctions 
in a closed form, which makes calculation of the propa- 
gator and determinant difficult. We thus once again wish 
to make an approximation that will enable us to calcu- 
late the propagato r and determinant. We note that the 
right hand side of (All) is approximately unity, allowing 
us to obtain 



k = 



2rnr 



Vn £ Z,n f 



(A12) 



for large values of k. The eigenfunctions ^>k, can then be 
approximated by a plane wave, and is given by 



(A13) 



with corresponding eigenvalue X n — (^j 2 -) +4|_E|. This 
implies, as is to be expected, that all the higher lying 
scattering states can be very well approximated by free 
particle states and that this only breaks down for the 
lowest lying scattering states, where the potential is im- 
portant. Note, however, that even the spectrum of the 



lower lying scattering states is well approximat ed by a 
free particle spectrum as the right hand side of ( |All| ) is 
approximately unity also in this case. The eigenfunctions 
are, however, distorted away from plane waves due to the 
presence of the potential. 

It is now possible to calculate the determinant using 
the above approximation so that 



!det(A± 1 )| = (3\E\) 



cx sinh 




4|£| 



(A14) 



where we have once agaia used the identity 1.143.1 from 
Gradshteyn and RyzhikcJ. Note that the quadratic term 
in ( A14| ) is due to the doubly degenerate bound state 
(from the symmetric and anti-symmetric eigenfunctions) , 
while the product over the continuum eigenvalues is 
raised to the fourth power since there is a fourfold degen- 
eracy in the continuum eigenfunctions (from the symmet- 
ric and anti-symmetric states, as well as the right moving 
and left moving plane waves). 

In calculating the propagator, we neglect the contri- 
bution of the bound state which gives only a small 
contribution to the propagator as it involves the prod- 
uct of two very well localized eigenfunctions evaluated a 
points which are well s epera ted. Taking into account only 
the scattering states, (A13), we find for the propagator 



A±(x,y) 



L 



E 



co S (^(x-y)) 



(A15) 
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